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We study the structural evolution of Sr 3 lr 2 07 as a function of pressure using x-ray diffraction. 

At a pressure of 54 GPa at room temperature, we observe a first-order structural phase transition, 
associated with a change from tetragonal to monoclinic symmetry, and accompanied by a 4% vol¬ 
ume collapse. Rietveld refinement of the high-pressure phase reveals a novel modihcation of the 
Ruddlesden-Popper structure, which adopts an altered stacking sequence of the perovskite bilayers. 

As the positions of the oxygen atoms could not be reliably refined from the data, we use density 
functional theory (local-density approximation-!-F+spin orbit) to optimize the crystal structure, 
and to elucidate the electronic and magnetic properties of Sr 3 lr 2 07 at high pressure. In the low- 
pressure tetragonal phase, we find that the in-plane rotation of the IrOe octahedra increases with 
pressure. The calculations further indicate that a bandwidth-driven insulator-metal transition oc¬ 
curs at ~20 GPa, along with a quenching of the magnetic moment. In the high-pressure monoclinic 
phase, structural optimization resulted in complex tilting and rotation of the oxygen octahedra, and 
strongly overlapping t 2 g and eg bands. The t 2 g bandwidth renders both the spin-orbit coupling and 
electronic correlations ineffectual in opening an electronic gap, resulting in a robust metallic state 
for the high-pressure phase of Sr 3 lr 207 . 


I. INTRODUCTION 

Understanding the metal-insulator transition (MIT) in 
transition metal oxides is pivotal in explaining the rich 
phenomenology that these compounds display [1]. Re¬ 
cently, strong spin-orbit coupling has been shown to lead 
to novel MITs and associated electronic and magnetic 
states in oxides containing a 5d element [2]. While the 
strongly correlated Sd oxides are described as Mott in¬ 
sulators, a range of mechanisms have been proposed to 
account for the insulating ground states of the 5d oxides 
[2-5]. In particular, it remains unclear how the metal- 
insulator transition evolves in the correlated, strong spin- 
orbit coupling regime relevant to these materials. 

In the Ruddlesden-Popper series of perovskite iridates 
Sr^+iRnOsn+i (where n is the number of Ir-0 layers, 
separated by Sr-0 layers), the large octahedral crystal- 
field splitting and spin-orbit coupling produce a Jeff = 
1/2 band for Ir^+ (5d^). For Sr 2 lr 04 (n = I), it was 
proposed that the narrow Jeff = 1/2 band is split by the 
modest electronic correlations to create a spin-orbit Mott 
insulator [3, 6, 7]. The ideal realization of the Jeff = 1/2 
state carries a magnetic moment of I/i^, however, ex¬ 
perimental values of 0.5/^^ [8, 9] suggest significant de¬ 
grees of itineracy. The perovskite iridates are therefore 
located in the intermediate-coupling regime, and the in¬ 
sulating state has been proposed to exhibit both Slater 
and Mott characteristics [5, 10, 11]. The bilayer com¬ 
pound Sr 3 lr 207 (n = 2) marginally retains the insulat¬ 


ing Jeff = 1/2 state [12-14]. The effective change in di¬ 
mensionality in the bilayer structure places Sr 3 lr 207 in 
close proximity to a metal-insulator transition, as Srlr 03 
(n = oo) is found to be metallic [4]. 

The charge gap of ^130 meV in Sr 3 lr 2 07 [15] can be 
closed by introducing mobile carriers through substitut¬ 
ing La^+ on Sr^+ sites [16-18]. This weakens electron 
correlations such that a conventional Fermi liquid state 
is recovered. In this context, application of pressure is 
an attractive control parameter, as it allows to continu¬ 
ously and isoelectronically tune the salient energy scales. 
Zhao et al. reported a second-order structural transi¬ 
tion in Sr 3 lr 207 at 14 GPa, accompanied by a reduction 
in resistivity [16, 19]. However, high-pressure transport 
measurements by Zocco et al. [20] have seen no evidence 
for this transition, and instead report a gradual reduc¬ 
tion of the electronic gap up to 104 GPa. While the 
authors do not report definitive evidence for a transition 
to a metallic state, they concede that pressure gradients 
could affect their data and that an insulator-metal tran¬ 
sition may occur at lower pressures. 

Here, we present a comprehensive study of the pres¬ 
sure dependence of the structure of Sr 3 lr 207 up to 61 
GPa, using powder x-ray diffraction (XRD) at ambient 
temperature and at 20 K. We identify a first-order struc¬ 
tural transition to a monoclinic symmetry accompanied 
by a volume collapse of 4% at 54 GPa at room tem¬ 
perature. Rietveld refinement suggests that the stack¬ 
ing sequence of the two perovskite bilayers is altered in 
the high-pressure phase, resulting in a novel modification 
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of the Ruddlesden-Popper structure. In contrast to the 
study by Zhao et al [19], no evidence for a second-order 
structural transition around 14 GPa was found either at 
room temperature or at 20 K. 

As the oxygen atoms scatter only weakly compared to 
the electron-heavy Sr and Ir atoms, a reliable refinement 
of the oxygen positions was not possible from the diffrac¬ 
tion data. We employed density functional theory [local- 
density approximation (LDA)+l/+spin-orbit (SO)] cal¬ 
culations to refine the structure of Sr 3 lr 207 at high pres¬ 
sure and to examine the evolution of the electronic and 
magnetic states. The calculations suggest that the in¬ 
plane rotation of IrOe octahedra increases with pressure 
in the tetragonal phase. The decrease of the Ir-0 bond 
distance, however, widens the t 2 g bands, such that an 
insulator-metal transition occurs at ^20 GPa. For the 
monoclinic high-pressure phase, optimization of the oxy¬ 
gen positions leads to complex rotations and tilting of 
octahedra, which considerably lowers the total energy 
with respect to the structure proposed by experiment. 
In this high-pressure structure, we find strongly over¬ 
lapping t 2 g and eg bands. Due to the wide t 2 g bands, 
both the spin-orbit coupling and electronic correlations 
become ineffectual in opening an electronic gap and we 
obtain a metallic state for the high-pressure phase. 


II. HIGH-PRESSURE X-RAY DIFFRACTION 
A. Experimental method 

We conducted high-pressure x-ray diffraction experi¬ 
ments at beamlines ID09a (European Synchrotron Radi¬ 
ation Facility) and 115 (Diamond Light Source). Sin¬ 
gle crystals of Sr 3 lr 207 were flux grown as described 
in Ref. 13 and characterized by SQUID magnetome- 
try. We used symmetrical diamond-anvil cells, fitted with 
Boehler-Almax diamonds with a culet diameter of 250 
fim. The gasket material was rhenium, indented to 30 
/im thickness, into which a 100-/im-diameter hole was 
laser-drilled. The Sr 3 lr 2 07 single crystals were ground 
into powders and loaded with Helium as the pressure 
transmitting medium. Pressure was measured in situ 
using the SrB 407 :Sm^+ fluorescence pressure scale [21]. 
The diffraction patterns were recorded on MAR345 and 
MAR555 detectors and integrated using the Fit2D soft¬ 
ware [22]. Full profile Rietveld refinements were per¬ 
formed using the FullProf software [23]. 


B. Results and discussion 

Diffraction patterns of Sr 3 lr 207 at high pressure are 
shown in Fig. 1. At ambient temperature in a he¬ 
lium pressure-transmitting medium, Sr 3 lr 207 remained 
in its tetragonal phase (space group M/mmm) up to 
pressures of 53 GPa. It has been argued that the ro¬ 
tations of the oxygen octahedra are correlated, resulting 



FIG. 1. (Color online). Integrated XRD spectra of Sr 3 lr 207 
at high pressure at ambient temperature, (a) Sr 3 lr 207 in its 
monoclinic phase (space group C2) at 61 GPa; (d) the corre¬ 
sponding structure, (b) Tetragonal phase of Sr 3 lr 2 07 (space 
group M/mmm) at 17 GPa; (c) the corresponding structure. 
In (a) and (b) the red circles are background-subtracted data 
points, the black solid line presents the Rietveld refinement, 
blue tick marks are the calculated Bragg reflections, and the 
residuals of the refinement are shown by the black line be¬ 
neath. Reflections from SrB 407 :Sm^^ are marked by an as¬ 
terisk. The inset in (a) shows the pressure evolution of small- 
angle reflections. In (c) and (d) the iridium atoms are shown 
in bronze, the oxygen atoms in red and the strontium atoms 
in purple. The structures in (c) and (d) were visualised via 
VESTA [24]. 


in a 90° twinned structure of pseudo-orthorhombic space 
groups Bhch and Acaa [25]. As our measurements were 
performed on powder samples, the tetragonal M/mmm 
structural model effectively averages over the two space 
groups. Rietveld refinements of our data confirm the 
M/mmm structural model to be adequate. Grucially, 
no second-order structural transition was observed near 
14 GPa. We used helium as the pressure-transmitting 
medium (at room temperature and 20 K), and not neon 
as in Ref. 19 (at 25 K). The discrepancy between our 
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studies may therefore arise from nonhydrostatic pressure 
conditions in the previous study. In our experiments, the 
pressure was determined by SrB 407 :Sm^+ fluorescence. 
No significant broadening of the SrB 407 :Sm^+ fluores¬ 
cence line shapes was observed, indicating that hydro¬ 
static pressure conditions were maintained. At pressures 
above 54 GPa, new reflections started to emerge. Fur¬ 
ther increase in pressure resulted in the complete disap¬ 
pearance of the /4/mmm reflections. The new reflec¬ 
tions were indexed in a C-centered monoclinic system, 
space group (72/m, (72, or Cm. Indexing the reflections 
in a higher symmetry proved not to be possible. This 
high-pressure phase remained stable up to 61 GPa, the 
maximum pressure reached in this study. Upon decom¬ 
pression, the monoclinic signature disappeared below 49 
GPa, and Sr 3 lr 207 reverted to its tetragonal phase. 

The high-pressure phase of Sr 3 lr 207 departs strongly 
from the Ruddlesden-Popper structure. Most notably, 
the reflection with the largest d spacing, (002) at 20 = 
2.4°, disappears [see inset Fig. 1(a)]. All new reflections 
possess smaller d spacings, indicating that the long c axis 
of tetragonal Sr 3 lr 207 is not retained. Figure 2 shows 
the refined lattice parameters as a function of pressure. 
In the tetragonal phase, the compression of Sr 3 lr 207 is 
anisotropic; the c axis exhibits a lower compressibility 
than the a=b axes. The c/a ratio thus increases with 
pressure by 3% to a critical value, where the system un¬ 
dergoes a structural phase transition. We indexed the 
new phase in a monoclinic symmetry, spacegroup (72/m, 
(72, or (7m, with lattice parameters a = 3.8351(4) A, 
b = 17.825(2) A, c = 3.5090(4) A, and f3 = 91.86(7)° 
at 61 GPa. This corresponds to a volume collapse of 
4% at the phase transition. A discontinuity of 2 A in 
the long c axis of the unit cell (labeled the b axis in the 
high-pressure phase) occurs at the structural transition. 
The lower compressibility along the c axis in the tetrag¬ 
onal phase appears to be alleviated by the phase transi¬ 
tion to monoclinic symmetry. A Vinet equation of state 
was fitted to the pressure-volume data of the tetrago¬ 
nal phase, yielding a bulk modulus of Bq = 157(4) GPa 
[Bq = 4.3(3)], in agreement with the literature [19]. 

To understand the transition at an atomic level, we at¬ 
tempted structural refinement of the high-pressure phase 
through Rietveld analysis. In general, the calculated in¬ 
tensities reproduce adequately the experimental intensi¬ 
ties when one of the perovskite bilayers of Sr 3 lr 207 is 
translated in plane, by half a unit cell [see Fig. 1(d)]. 
This may explain the discontinuity in the lattice parame¬ 
ter of the long axis of the unit cell. While in its tetragonal 
phase the compressibility along the c axis appears to be 
limited, this is alleviated in the monoclinic phase by an 
alteration in the stacking sequence of the perovskite bi¬ 
layers. Furthermore, the Sr^+ atoms are arranged within 
a single plane, parallel to the a axis in the unit cell. To 
accommodate this configuration, the lattice parameter a 
increases, such that the in-plane lattice constants (a, c) 
of the monoclinic structure differ by 10%, and the angle 
spanned between a and c is 91.8°. 


(a) Pressure [GPa] (b) Pressure [GPa] 



Pressure [GPa] 

FIG. 2. (Color online). Equation of state of 81311207 . (a) 
and (b) Lattice parameters as a function of pressure at 298 
K. The inset in (a) shows the normalized lattice constants 
as a function of pressure at 20 K. (c) Equation of state of 
Sr 3 lr 2 07 at 298 K. The solid line is a Vinet-type fit to the 
observed equation of state of Sr 3 lr 207 . The inset in (c) shows 
the c/a ratio as a function of pressure at 298 K. 


This high-pressure structure could be described in the 
space group (72/m. The presence of a mirror plane per¬ 
pendicular to the b axis minimizes the number of free 
parameters to two y coordinates for Ir and Sr (the sec¬ 
ond Sr atom is fixed at ^ = 0.5). However, the fit of 
the calculated diffraction pattern was deemed unsatisfac¬ 
tory. We therefore decided to remove the mirror plane 
and adopt (72 symmetry. This allowed to move atoms 
independently along the y axis, through two adjustable 
y coordinates for both Ir and Sr atoms (here, the third 
Sr atom remained at ^ = 0.5). This provided the best fit 
to the data. The weighted profile Rietveld R factor 
of the monoclinic structure at 61 GPa then became com¬ 
parable to the one obtained from the Rietveld refinement 
of the tetragonal phase at 17 GPa. Thus our proposed 
structural model for the high-pressure phase provides as 
good a description of the data as the /4/mmm struc¬ 
ture at lower pressures. We note that the visually worse 
fit of the calculated intensities at 61 GPa at least partly 
originates from the difficulty of accurately describing the 
background at very high pressures. 

As the oxygen atoms scatter only weakly compared 
to the electron-heavy Ir and Sr atoms, a precise refine¬ 
ment of the oxygen positions was not possible for the 
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high-pressure phase. Small deviations from the proposed 
positions [as shown in Fig. 1(d)] through rotations or tilt¬ 
ing cannot be excluded or confirmed within the quality 
of our data. 


III. RESULTS OF THE FIRST-PRINCIPLES 
CALCULATIONS 

A. Numerical methods 


To investigate the pressure dependence of the elec¬ 
tronic and magnetic properties of Sr 3 lr 207 we performed 
first-principles calculations with two complementary ap¬ 
proaches: the full-potential linearized-augmented-plane 
wave (FP-LAPW) method, as realized in the Elk pack¬ 
age [26]; and the plane wave pseudopotential method, 
as implemented in the QUANTUM-ESPRESSO simu¬ 
lation package [27]. In these calculations we used the 
LDA, and its popular refinement taking into account 
the on-site Coulomb interaction and spin-orbit coupling 
(LDA+U+SO) [28]. In pseudopotential calculations, we 
set an energy cutoff of 400 eV in the plane-wave ba¬ 
sis construction and the energy convergence criteria of 
10“^ eV. Eor the Brillouin-zone integration, a 6 x 6 x 6 
Monkhorst-Pack mesh was used. 
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FIG. 3. (Color online). (Left) Structure of Sr 3 lr 207 in the 
ambient, tetragonal phase (space group lA/mmm) with oxy¬ 
gen disorder. (Middle) Bbcb structure of Sr 3 lr 207 used in 
first-principles calculations. The difference between blue and 
gray IrOe octahedra is the octahedral rotation around the z 
axis. The red and purple spheres denote oxygen and stron¬ 
tium atoms, respectively. (Right) xy-plsiue projection of the 
Bbcb structure. Crystallographic plots were made with the 
VESTA software [24]. 


B. Structural and electronic properties at ambient 
pressure 

We first discuss the LDA and LDA-fSO results for the 
ambient pressure phase of Sr 3 lr 207 . In these calcula¬ 
tions, we started from the /4/mmm crystal structure 
with oxygen disorder [29] and then doubled the planar 
basis lattice vectors and deleted specific planar oxygen 



FIG. 4. (Color online). The partial densities of states calcu¬ 
lated for Sr 3 lr 207 in the Bbcb structure at ambient pressure 
using the LDA method. 


atoms in the supercell, to achieve the Bbcb space group 
(Fig. 3). 

From the calculated partial densities of states pre¬ 
sented in Fig. 4, the splitting between t 2 g and Cg states 
can be estimated to be about 3.3 eV, in good agreement 
with the experimental value of 3.5 eV [30]. We find that 
the bd Ir states of xy character are filled and the yz (xz) 
orbitals are partially occupied. The Cg band is fully lo¬ 
cated in the unoccupied part of the electronic spectrum 
due to much stronger hybridization with oxygen states. 
The obtained level picture agrees with symmetry anal¬ 
ysis of the Sr 3 lr 207 system [8]. As we show below, the 
degeneracy of xz and yz states plays an important role 
for the orbital magnetism. 

The comparison of the LDA and LDA+SO band struc¬ 
tures (Fig. 5, left) shows that the spin-orbit coupling is 
responsible for the formation of four doubly degenerate 
nearly separated narrow bands at the Fermi level. The 
corresponding bandwidth can be estimated as 0.5 eV. 
Given the large crystal field, an effective total angular 
momentum approach can be used to construct a mini¬ 
mal microscopic model of the insulating ground state of 
Sr 3 lr 207 , as proposed in Ref. 3. However, in contrast to 
the ideal Jeff = 1/2 state described in Ref. 3, we have 
a more complicated situation. The splitting of the t 2 g 
manifold results in a mixing of Jeff =1/2 and Jeff = 3/2 
bands, which should be taken into account within a mi- 
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FIG. 5. (Color online). The band structures of Sr 3 lr 207 at ambient pressure obtained from LDA (left), LDA+SO (middle) 
and LDA+f/+SO (right) calculations using the Bbcb structure. The coordinates of the symmetry points, G=( 0, 0, 0), Z=( 0, 
|, 0), T=( 0, |, ), Y=( 0, 0, -|^), S=( |, 0, ^), R=( |, |, are in units of where a = 5.5 A and c = 20.87 A. 


croscopic analysis [31]. 

We employed the LDA+7/+SO method to first repro¬ 
duce the insulating ground state and magnetic proper¬ 
ties of Sr 3 lr 207 at ambient pressure. The obtained band 
structure is illustrated in Fig. 5. In our simulations, 
the on-site Coulomb interaction U and intra-atomic ex¬ 
change parameter Jh were chosen to be 1.8 and 0.3 eV, 
respectively, as discussed in Ref. 32. The electronic gap 
obtained from the calculations is 110 meV, in reason¬ 
able agreement with the experimental estimate of 130 
meV [15]. The calculated spin and orbital magnetic mo¬ 
ments are = 0.14//^ and Ml = 0.24//^, respectively, 
in good agreement with the total magnetic moment of 
0.52(8)/iB measured by neutron diffraction [9]. It is in¬ 
teresting to note that the calculated orbital moment of 
Sr 3 lr 207 is comparable to the one calculated for Sr 2 lr 04 
[33], while the spin moment is two times larger. 

We have simulated different c-oriented magnetic con¬ 
figurations of Sr 3 lr 207 : ferromagnetic, and antiferromag¬ 
netic of G, G, and A types. We find that the c-axis 
antiferromagnetic configuration of G-type minimizes the 
total energy of the system, in agreement with experi¬ 
ments [12, 13]. The calculations for other types of the 
magnetic order revealed strong changes in electronic and 
magnetic properties of Sr 3 lr 207 . In general, the mag¬ 
netic moments and energy gap values are suppressed in 
comparison with those obtained for G-type structure. 

In turn, the in-plane antiferromagnetic solution with 
weak ferromagnetism also corresponds to the excited 
states of the Sr 3 lr 207 system. The corresponding en¬ 
ergy difference between in-plane and out-of-plane G-type 
solutions is 353 meV (per unit cell containing four Ir 
atoms). The calculated orbital and spin magnetic mo¬ 
ments of Ir atom in the in-plane magnetic configuration 
are 0.186//^ and 0.063//^, respectively. The canting an¬ 
gles of the orbital and spin moments are equal to 15.2° 
and 18.4°, respectively. As in the case of Sr 2 lr 04 [33] 
we found a strong contribution from the oxygen atoms 
to the resulting net magnetic moment. 


C. Structural and electronic properties at high 
pressure 

To understand the evolution of the electronic and mag¬ 
netic properties of Sr 3 lr 207 as a function of pressure, we 
have performed LDA+7/+SO calculations for the Bbcb 
structure, using the experimental lattice parameters (Fig. 
2). For each pressure point, the crystal structure was op¬ 
timized. The resulting structural, electronic, and mag¬ 
netic properties of Sr 3 lr 207 as a function of pressure are 
shown in Table 1. 

Our calculations reveal that the compression of the 
unit cell is accompanied by an increase of the in-plane 
rotation angle of the IrOe octahedra, while the Ir-0 dis¬ 
tance decreases. Intuitively, this can be understood from 
the experimentally observed anisotropic compression of 
Sr 3 lr 207 , where the compressibility along the a and b 
axes is higher than that along the c axis. At ambient 
pressure, the rotation angle is about 11.8°, while the 
high-pressure solution at 50 GPa reveals an octahedral 
rotation of 16.3°. 

For Sr 2 lr 04 and Ba 2 lr 04 , it was demonstrated that 
the in-plane rotation of the IrOe octahedra with fixed Ir- 
O distance leads to a more localized electronic structure 
[33]. Therefore the width of the t 2 g band was found to 
be narrower for Sr 2 lr 04 than for Ba 2 lr 04 . For Sr 3 lr 207 , 
applying pressure also decreases the Ir-0 distance, which 
leads to a strong increase of the LDA bandwidth for t 2 g 
states. Figure 6 shows that the spectral function for de¬ 
generate xz and yz states retains its structure at high 
pressure; however, the width of the spectral function no¬ 
tably increases at 50 GPa. 

In our calculations, the insulating state of Sr 3 lr 207 
is therefore very sensitive to the value of the pressure. 
Table I shows the obtained dependence of the electronic 
gap on pressure. Above 20 GPa, we find a closure of the 
electronic gap and a transition to a metallic state. Ex¬ 
perimentally, Li et al. reported a transition to a nearly 
metallic state at 13.2 GPa [16] with a vanishingly small 
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TABLE L Pressure dependence of the structural, electronic, 
and magnetic properties of Sr 3 lr 207 , using the Bbcb struc¬ 
ture. The results were obtained with LDA+f/+SO calcula¬ 
tions, as described in the text, a is the angle of the in-plane 
Ir-O-Ir bond. E^ap is the value of the energy gap. The lat¬ 
tice parameters corresponding to the pressure of 60 GPa were 
obtained by extrapolation of the experimental data. 


P (GPa) 

a 

dir-o (A) Egap (meV) Ms (ys) 

Ml (ys) 

0 

156.4 

1.99 

no 

0.14 

0.24 

10 

155.4 

1.95 

62 

0.14 

0.21 

20 

153.1 

1.93 

- 

0.11 

0.14 

30 

150.9 

1.91 

- 

0.07 

0.08 

40 

148.4 

1.89 

- 

0.05 

0.05 

50 

147.4 

1.88 

- 

0.04 

0.04 

60 

144.6 

1.87 

- 

0.04 

0.03 



EIG. 7. (Golor online). Total densities of states of Sr 3 lr 207 
in the Bbcb structure, calculated at different pressures using 
the LDA+f/+SO approach. Solid, dashed and dotted lines 
correspond to ambient pressure, P = 20 GPa and P = 30 
GPa, respectively. 


electronic gap. However, high-pressure transport mea¬ 
surements by Zocco et al [20] suggest that an electronic 
gap remains up to at least 30 GPa. The calculated den¬ 
sities of states of Sr 3 lr 207 near the Fermi level as a func¬ 
tion of pressure are shown in Fig. 7. There is a pseudo¬ 
gap state at P = 20 GPa; for larger pressures we obtain 
a metallic solution with non-zero density at the Fermi 
level. 

Another interesting effect evident in our calculations 
is a suppression of the magnetic moment at high pres¬ 
sure. Analysis of the calculated occupation matrices 
shows that the strongest contribution to the orbital mag¬ 
netic moment of Ir comes from the 5d states with mag¬ 
netic quantum number m/ = ± 1 that correspond to 
partially-filled, degenerate xz and yz states in the LDA 
spectrum. This contribution is suppressed under pres¬ 
sure due to strong hybridization with oxygen states. An 
indication of the hybridization amplification is the LDA 
spectral function spreading at high pressures (Fig. 6). 



FIG. 6. (Golor online). Gomparison of the partial densities 
of states of Sr 3 lr 2 07 in the Bbcb structure, calculated using 
the LDA method at ambient pressure (blue solid line) and 50 
GPa (red dotted line). 


D. High-pressure phase of Sr 3 lr 2 07 


Analysis of the experimental data established a struc¬ 
tural transition to a monoclinic phase at 54 GPa. Since 
a precise refinement of the oxygen positions from the ex¬ 
perimental data was not possible, we explored different 
arrangements of the oxygen octahedra to optimize the 
crystal structure of Sr 3 lr 2 07 in the C2 phase. We sim¬ 
ulated supercells with and without rotations of the oxy¬ 
gen octahedra (Fig. 8). To examine the nature of the 
structural transition, we also compare our results with 
a hypothetical Bbcb-type structure at 61 GPa, obtained 
through an extrapolation of the lattice parameters. 

We first discuss the dependence of the refined struc¬ 
ture without oxygen rotations on the choice of the first- 
principles method. The structures obtained with the 
LDA, LDA+SO, and LDA+1/+SO approaches were com- 



FIG. 8. (Golor online). High-pressure structures of Sr 3 lr 207 
used in this work. (Left) The structure without rotations of 
the oxygen octahedra, as determined by experiment. (Right) 
The optimized structure corresponding to the ground state 
at high pressure, P = 61 GPa. The octahedra belonging to 
the A (D) and B (G) layers rotate around out-of-plane and 
in-plane axes, respectively. 
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TABLE IL Comparison of Ir-0 distances (in A) of the exper¬ 
imental and refined high-pressure C2 structure of Sr 3 lr 207 , 
without rotations of the oxygen octahedra (Fig. 8 , left), at 61 
GPa. The calculations were performed with the experimental 
lattice parameters, using different first-principles methods. 


bond 

expt. 

LDA 

LDA+SO LDA+H+SO 

Ir-Ol (in-plane) 

1.7545 

1.7589 

1.7597 

1.7562 

Ir-02 (in-plane) 

1.7545 

1.7589 

1.7597 

1.7562 

Ir-03 (in-plane) 

1.8449 

1.8806 

1.8811 

1.8804 

Ir-04 (in-plane) 

1.8449 

1.8937 

1.8974 

1.8989 

Ir-O5 (out-of-plane) 

1.9175 

1.9223 

1.9222 

1.9189 

Ir-06 (out-of-plane) 

1.9175 

1.9223 

1.9222 

1.9189 


pared (Table II). In all cases, we used the experimen¬ 
tal lattice parameters. A tendency to form nonequiv¬ 
alent in-plane bonds is evident. Importantly, the dif¬ 
ference in interatomic distances obtained with LDA, 
LDA+SO and LDA+U+SO methods is between 10“^ A 
and 10“^ A. Such a weak dependence of the optimized 
structure on the choice of the first-principles method 
was found with both the full-potential Elk package and 
the plane-wave pseudopotential approach, realized in the 
Quantum-Espresso simulation package. The total densi¬ 
ties of states shown in Fig. 9 corresponds to a metallic 
state. 

The LDA solution for the C2 structure without rota¬ 
tions of the oxygen octahedra (as determined by experi¬ 
ment; Fig. 8, left) is higher in energy by about 3 eV (per 
unit cell containing two Ir atoms) than the Bbcb struc¬ 
ture at 61 GPa. This largely originates from the more 
compressed volume of the unit cell of the C2 structure 
(240 A^), compared to the (extrapolated) volume of the 
Bbcb structure at 61 GPa (250 A^). 

We now consider the optimized C2 structure of 
Sr 3 lr 207 , which adopts complex rotations of the oxy¬ 
gen octahedra (Fig. 8, right). In this case, the LDA 
calculations revealed that the lowest total energy of the 
system corresponds to a geometry in which the octahe¬ 
dra from the A and D layers rotate in plane (around the 



FIG. 9. (Color online). Comparison of the total densities 
of states of the refined high-pressure structure of Sr 3 lr 2 07 
without rotation of the oxygen octahedra, obtained from 
LDA (black solid line), LDA+SO (blue dotted line) and 
LDA+I/+SO (red dashed line) calculations. The Fermi level 
corresponds to zero energy. 




FIG. 10. (Color online). Top (left panel) and side (right 
panel) views of the oxygen octahedra belonging to the B(C) 
plane of Sr 3 lr 207 in the optimized high-pressure structure, as 
shown in the right panel of Fig. 8 . 


long b axis) by about 9°. This in-plane rotation is simi¬ 
lar to the one observed in the ambient pressure structure. 
Conversely, the rotation axis for the octahedra from the 
B and C planes is out of plane (tilting), as shown in Fig. 
10. We obtained a tilting angle of about 20°. We also 
found a tendency for forming short (1.85 A) and long 
(1.94 A) bonds within the IrOe octahedra (Table III). 

This optimized high-pressure structure of Sr 3 lr 207 
(Fig. 8, right) is 1.4 eV lower in energy (per unit cell 
containing two Ir atoms) than the experimentally deter¬ 
mined high-pressure structure without rotations. We be¬ 
lieve that further lowering of the total energy is possible 
by rotating the IrOe octahedra in the A and D planes 
around an in-plane axis. However, this requires an in¬ 
crease in the size of the supercell and optimization of the 
geometry. 

It is important to analyze the difference in the 
electronic structure of the high-pressure structure of 
Sr 3 lr 207 . Figure 11 compares the total densities of states 
of the compressed Bbcb and optimized C2 structures, ob¬ 
tained with LDA. The plane-dependent rotation of the 
oxygen octahedra and compression of the unit cell in the 
C2 structure lead to a complete delocalization of the 5d 
Ir states in Sr 3 lr 207 . No electronic gap between the t 2 g 
and Cg bands remains. 

Our LDA+SO and LDA+U+SO calculations for an 
optimized C2 structure with a complex rotation of the 
oxygen octahedra have demonstrated that neither spin- 
orbit coupling nor on-site Coulomb interaction can signif¬ 
icantly affect the electronic structure. The corresponding 
densities of states are presented in Fig. 12. The system 
remains metallic, which still follows from the increased 
t 2 g bandwidth. 


TABLE III. Ir-O distances of 81311207 in the optimized high- 
pressure structure with rotations of the oxygen octahedra, as 
shown in Fig. 10. 

plane d/r-o (in-plane) d/r-o (apical) 

A 1.85 - 1.94 A 1.89 - 1.91 A 

B 1.87 - 1.92 A 1.91 - 1.91 A 
















FIG. 11. (Color online). LDA total densities of states of 
Sr 3 lr 207 in the compressed Bbcb structure (blue dotted line) 
and optimized C2 structure (red solid line) at P = 61 GPa. 
The Fermi level corresponds to zero energy. 



FIG. 12. (Color online). Comparison of the total densities 
of states of the optimized C2 structure of Sr 3 lr 207 obtained 
from LDA (black solid line), LDA+SO (blue dotted line) and 
LDA+P+SO (red dashed line) calculations . The Fermi level 
corresponds to zero energy. 

IV. CONCLUSION 

We have studied the structural evolution of Sr 3 lr 207 
with powder x-ray diffraction to unprecedented pres¬ 
sures. We observe a novel first-order structural transi¬ 
tion at 54 GPa, preceded by a strongly anisotropic com¬ 
pression, characterized by an increasing c/a ratio with 
pressure. The high-pressure phase could be indexed in a 
monoclinic symmetry (space group (72/m, (72, or Cm). 
Rietveld refinement of the diffraction data suggests that 
the stacking sequence of the perovskite bilayers is altered 
in the high-pressure structure. A precise refinement of 
the oxygen positions was, however, not possible due to 


the relatively weak x-ray scattering signal compared to 
Ir and Sr atoms. 

We used first-principles LDA+(7+SO calculations to 
both optimize the structure and examine the electronic 
properties of Sr 3 lr 207 as a function of pressure. Below 
the high-pressure structural transition, we find that the 
anisotropic compression observed by experiment is ac¬ 
companied by an increase in the in-plane IrOe octahedral 
rotation. The simultaneous decrease of the Ir-0 bond 
length however increases the t 2 g bandwidth and gradu¬ 
ally closes the band gap. In our calculations, we find 
an insulator-metal transition at ^20 GPa. This value is 
comparable to the transition to a nearly metallic state 
at 13.2 GPa observed by Li et al [16], but contrasts 
with measurements by Zocco et al. [20], which show an 
electronic gap remaining up to at least 30 GPa. 

For the high-pressure phase, structural optimization 
of the oxygen position resulted in complex in-plane 
and out-of-plane rotations of the octahedra. This 
considerably lowered the total energy with respect to the 
structure obtained from experiment. In the optimized 
high-pressure structure, the t 2 g and Cg bands strongly 
overlap. The shortened Ir-0 bond length does not allow 
an electronic band gap to open, even in the presence 
of strong spin-orbit coupling and on-site Coulomb 
interactions. This suggests a robust metallic state for 
the high-pressure monoclinic phase of Sr 3 lr 2 07 . 
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